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Abstract 

Approximate solutions for potential flow past an axisymmetric slender body 
and past a thin airfoil are calculated using a uniform perturbation method and 
then compared with either the exact analytical solution or the solution 
obtained using a purely numerical method. The perturbation method is based 
upon a representation of the disturbance flow as the superposition of 
singularities distributed entirely within the body, while the numerical 
(panel) method is based upon a distribution of singularities on the surface of 
the body. It is found that the perturbation method provides very good results 
for small values of the slenderness ratio and for small angles of attack. 
Moreover, for comparable accuracy, the perturbation method is simpler to 
implement, requires less computer memory, and generally uses less computation 
time than the panel method. In particular, the uniform perturbation method 
yields good resolution near the regions of the leading and trailing edges 
where other methods fail or require special attention. 
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Introduction 


Many different methods have been developed to obtain an approximate 
solution to potential flow problems involving fixed finite bodies. For a few 
very special cases, analytical methods can be used to obtain an exact solution 
to the exterior flow problem. However, in most cases, usually due to the 
geometric complexity of the problem, approximate methods must be employed. 
Two of the more widely used approximate methods are perturbation methods and 
numerical methods. Perturbation methods take advantage of the fact that many 
shapes of practical interest are either slender or thin and hence certain 
simplifications can be made in the problem formulation. Many numerical 
methods are based on the idea of approximating the body as a collection of 
"simpler" shapes (e.g., panels) and then satisfying an appropriate boundary 
condition on each of these simpler shapes. The primary purpose of this paper 
concerns the results of comparing the approximate solutions obtained by a 
uniform perturbation method and a standard numerical method (and exact 
solutions, where they exist) to a variety of potential flow problems involving 
either a slender body of revolution or a thin airfoil. 

The perturbation method we shall use is based upon an approach suggested 
by Handelsman and Keller [1] and developed more fully by Geer and Keller [2] 
and Geer [3,4]. This method represents the part of the potential due to the 
body as the superposition of potentials due to distributed point 
singularities. These singularities lie entirely within the body. The 
boundary condition on the surface of the body then leads to an integral 
equation for the density of the singularity distribution. This equation can 
be solved asymptotically as the slenderness ratio of the body approaches 
zero. By demanding that the resulting asymptotic expansion be regular, the 
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extent of the singularity distribution within the body can also be 
determined. Thus, this method leads to an approximate solution for the flow 
field which is uniformly valid over the entire surface of the body, including 
the ends (e.g. , nose and tail) of the body. Hence, as we shall see, it can 
easily be used to compute the surface velocity (and hence the pressure 
coefficient) on the body surface. 

The numerical method we shall use is commonly referred to as a "panel 
method (see, e.g., Maskew [5]). Here the body is approximated by a finite 
number of panels (or line segments, in two dimensions). Then the flow due to 
the presence of the body is represented by a distribution of singularities 
over these elements, with a polynomial variation of the density of the 
distribution over each element. The boundary condition on the body surface 
then leads to a system of algebraic equations for the polynomial coefficients. 

In Sections 2 and 3, below, the specific problems that we shall consider 
are formulated, and the uniform asymptotic solutions are presented. Surface 
pressure coefficients computed from the perturbation and numerical methods are 
presented in Section 4, while our results are discussed in Section 5. 


2. Flow Past an Axisymmetric Slender Body 

In this section and the next, we shall formulate two classes of potential 
flow problems that we shall consider in detail. In particular, these include 
flow past a slender body of revolution and the flow over a thin two- 
dimensional airfoil. For each class of problems, we will summarize the 
general uniform perturbation method of solution and present the relevant 
formulas. We shall then apply the results to some specific examples. 
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Consider the problem of determining the steady potential flow past a 
slender body of revolution (see Figure 1). We fix a cylindrical coordinate 
system (r,0,x) within the body, with the x-axis coinciding with the axis of 
the body and the origin at one end of the body, and assume, for simplicity, 

that the incident flow is a uniform flow with speed U parallel to the axis 

of the body. Then the surface of the body can be described by r = e/S(x) 
for 0 < x <_ 1. Here e is called the slenderness ratio of the body and is 

defined as the ratio of the maximum radius of the body to its length L. The 

function S(x), which satisfies max S(x) = 1 for 0 £ x < 1, is a prescribed 
function which describes the shape of the body, with ire S(x) being the 
cross - sectional area of the body at the station x. Here we have non - 
dimensionalized all lengths by referring them to L. 

To determine the flow about this body, we represent the part of the 
potential due to the body as the superposition of potentials due to point 
sources distributed along the axis of the body and lying entirely inside the 
body. Thus, we write the potential <|> as 


4 = U (x - 


f(5,e) 


4ir 


d£ 


“ J (X - O 2 * r 2 


( 2 . 1 ) 


Here f(5,e) is the (unknown) density of the source distribution, while a 
and B, which depend upon e and are also unknown, determine the extent of 
the distribution within the body. They satisfy the inequalities 
0 < a < B < 1; (see Figure 1). 

The potential defined by (2.1) satisfies Laplace's equation outside the 
body and reduces to the uniform flow potential at infinity. The boundary 
condition that the normal component of the flow velocity must vanish on the 
surface of the body, when used with (2.1), leads to the requirement 
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e 2 ^(x) - / — (x - — f(g.e)d£. (2.2) 

“i(x-0 2 + e 2 S(x) 

Equation (2.2) is an integral equation for the density of the source 
distribution, from which f(S,e), as well as a and 0, can be determined. 

For small values of e, the integral operator on the right side of (2.2) 
can be expanded asymptotically as a power series in e 2 and a power series in 

e multiplied by log(e 2 ). The coefficients in these series can be expressed 

as certain linear operators applied to the density function f(5,e); (see, 
e.g., Handelsman and Keller [1] for details). The form of this expansion 
suggests that the solution for f(£,e) can be expressed in the form 

f(x,e) = e f 1Q (x) + e 4 f 2Q (x) + e 4 log e 2 f 2J (x) + e 6 f 3Q (x) 

+ e 6 log e 2 f 31 (x) + e 6 (log e 2 ) 2 f 32 (x) + 0(e 8 log e 2 ), 

(2.3) 

where each of the functions f n>m (x) is independent of e . By substituting 
(2.3) into the expanded form of (2.2) and equating the coefficients of like 
terms in e> the functions f n>m (x) can be determined recursively. In 


particular, 


f 10 (x) = tfS'(x ) , 


£ / \ tt d ( ss ss cc // _l co"i ^ ^0 i 

f 20 (x) ’ '?S — -T^T- SS + SS l0g l STI7— J 


1-x 


-2 


+ S / (S'(x + v) - S'(x) - vS"(x)}v dv 

0 


- S / {S'(x - v) - S'(x) + vS"(x)}v 2 dv j 

0 


f 21 (x) = i f 30 (x) = 0 


f 31 ( x) = 


dx (Sf 20 ) 


1 fL 

8 dx 


( g2 £ ;a') - th w ** 1 / — 1 - s i 


dx v 21 


+ f 2l(| " + S J {f 21 (x " v) " f 21 (x) + vf 2i (x)}v ~ 2 dv 

1— X 

- S / { f 2j(x + v) - f 21 (x) “ vf 21 ( x ^ v 2 dv 


f 32 (x) = l£d^ < S 


j2 

^—7 (SS")}. 
dx 


(2.4) 


By demanding that each of the coefficient functions f n>ra (x) is regular 
on 0 _< x _< 1, Handelsman and Keller [1] were able to show that a(e) and 
3(e) can be expressed in the form 



a 


a e 
n 


(2.5) 


a(e) = l 
n=l 


2n 


0 = 0(e) 


1 ’ 6 n e 
n=l 


2n 


where the coefficients and 8 n can be expressed in terms of S(x). In 

particular, they showed that 


a 


C 1 2 
4 e 


C 1 c 2 _4 , c l^ c l c 3 + 2 c 2^ 6 

16 e + 64 e 


+ 0(e 8 ) , 


c . = 


S (j) (0)/j! 


( 2 . 6 ) 


d l 2 d l d 2 4 
8 " 1 ~ — e + 16 - G 


! l( d l 


64 


d 3 +2d 2) 6 +n ,8, 

e + 0(e ), 


d. 

J 


( _l)j S (j) (l)/j! 

(2.7) 


Thus, once S(x) is specified, a and 8 can be determined from (2.6)- 
(2.7) and then used in (2.1), along with the expression (2.3) for f(x,e), to 
yield the asymptotic expansion of the potential with an accuracy of 
0(e 8 log e 2 ) . 

As applications of these results, we consider two specific examples. 
These are: 

1) Ellipsoid of revolution, for which 


S(x) = 4x(l - x), (Figures 2a, b); 


2) "Dumbbell" shaped body, for which 


S(x) = 4bx( 1 - x)[l - bx( 1 - x)], b > 2, (Figures 3a, b). 
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When these specific formulas for S(x) are used in (2.4), (2.6), and (2.7), 
and then these results are used in (2.1), an explicit expression for the flow 
potential is obtained for each example. In Section 4 we shall use these 
formulas to compute the pressure coefficient on the surface of each body and 
then compare these results with the pressure coefficient found by using an 
exact solution and/or the panel method. 


3. Flow Past a Thin Airfoil 

Consider now the problem of determining the steady potential flow about a 
thin airfoil which is immersed in an otherwise uniform flow field. We 

introduce an x,y-cartesian coordinate system fixed in the airfoil (see Figure 
4) and let the surface of the airfoil be described by y = e[c(x) ± /S(x) J 

for 0 < x < 1. (Here we have again nondimensionalized x and y by 

referring all lengths to the chord L of the airfoil.) The parameter e is 
called the slenderness parameter of the body, where C(x) and S(x) are 

prescribed functions describing the camber line and thickness profile, 
respectively, of the airfoil. In particular, for an airfoil that is symmetric 
about the x-axis, C(x) = 0. 

To determine the flow potential <f> about the airfoil, it is convenient to 
think of <{> as the real part of a complex potential 4>(z), which is an 
analytic function (outside the airfoil) of the complex variable z = x + iy. 
Then we represent the part of the potential due to the body as the 
superposition of potentials due to (complex) point sources distributed along 
an arc which lies entirely within the body (see Figure 4). Thus, we set 


( 3 . 1 ) 


£ 

*(z) = u{e" iY z “ 27 / log(z - 5)f(C,e)d5}, 

a 

where f(5,e) is the (unknown) density of the source distribution while a 
and 3, which are also unknown, determine the extent of the distribution 
within the body. The parameter y represents the angle of the uniform stream 
makes with the positive x-axis (i.e., the angle of attack). The real part of 
this expression satisfies Laplace's equation outside the airfoil and reduces 
to the uniform flow potential at infinity. 

The only remaining condition to be satisfied is the vanishing of the 
normal component of the velocity on the body's surface. This condition, when 
used with ( 3 . 1 ), leads to a pair of coupled integral equations for the 
(complex) source density f(5,e), from which a and 3 can also be 

determined. The case of a symmetric airfoil (C(x) = 0 ) was analyzed by Geer 
and Keller [ 2 ], while the general case of a cambered airfoil was studied in 
detail by Geer [ 3 ]. Using the results of these investigators, we find that 
for small values of the slenderness parameter e, 

f ( z » £ ) 1 {f n (2) + f,(z)e + f (z)e 2 + 0(e 3 ) }, (3.2) 

AS - z)(z - a) U 1 2 
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where , 
f 0 (x) 

Ref j (x) 
Imf j(x) 

f j(x) 

Ref 2 (x) 
Imf 2 (x) 
f 2 (x) 


if— + (2x - l)siny) 

TT 


^s(x) ~ {ZitS (x)cosy + ReL^Cf^^x))} 

* 4 - f f ? i <si - 2 9W9 

{ ImL^(f q(x) ) - 4irC'(x)cosy} 

1 2 a n C=sin 2 9 

_ x j 2 cos e ._ {imL^f (?)) - 4itC'(5)cosy} d9 

it 0 sin 0 - x 5=x 

- 27 -J X( S(x) ~ t ReL l‘ £ l(*» + R ^2 (£ 0 (l!)> ) 

r - 

f„(x) + — - - J 2 {f „(sin 2 0) + (S'(O)cos 2 0 + S'(l)sin 2 0)}d0 

2 7T TT ' Z / 


2 ^- {lmL*j*(f j(x)) + ImL 2 (f Q (x))} 


i 2 0 _ ?=sin 0 

_ x j 2 cos 9 ( ImL 0 (f (5)) + i mL 2(f (J-))} d0 (3.3) 

* Z 0 sin 0 - x 1 1 * 5=x 


and a and 3 are given by 
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° = a(e) e 2 -| Cl ^ e 3 -J. d 1 (d 2 + 4c^)e 4 + 0(e 5 ) 

B = 3(e) = 1 “ Tj- e 2 - dj e 3 + yi d^d^ + 4c 2 )e 4 + 0(e 5 ) 

Cj = C C3) (0)/j ! ; d. = S (j) (0)/j!; Z. = (-l) j C (j) (l)/j!; d ' = (-l) j S (J) (l)/j ! 

(3.4) 

In (3.3), Lj are certain linear operators which are defined by Geer [3]. The 
constants I\ are related to the total circulation T(e) about the airfoil 
by f(e) - Tq + e + r 2 e + 0(e ). If the airfoil has a sharp trailing 
edge at z = 1, then Geer [3] has shown that the Kutta condition leads to the 
requirement that Imf-j(l) = 0 for each j. Hence the constants r are 

J j 

determined uniquely for this case from the formulas (3.3). Thus, once C(x) 
and S(x) have been specified, they can be used in (3.2)-(3.4) and then these 
results inserted into (3.1) to yield an approximate expression for the flow 
potential which has an accuracy of O(e^). 

As applications of these results, we consider two specific examples. 
These are the: 

1) Symmetric second-order Joukowski airfoil, for which 

C(x) = 0 and S(x) = x(l - x)^; (Figures 5a, b). 

2) Cambered second-order Joukowski airfoil, for which 


C(x) = ax(l - x) and S(x) = x(l - x)^; (Figures 6a-d). 
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In the next section, we will use the formulas above to compute the surface 
pressure coefficient for these examples, with different angles of attack. 


4. Comparison of Perturbation and Numerical Results 

In this section we shall compare the results of using the perturbation 
formulas of the previous two sections and a low-order panel method code (see, 
e.g., Maskew [5]) to compute the pressure coefficient C p on the surface of a 
slender body of revolution or a thin airfoil. These results will also be 
compared with those obtained from an exact analytical solution to the flow 
problem, when such a solution is available. In particular, using Bernoulli's 
equation, we find 

C p = (2/ pU 2 )(p - pj = 1 - U -2 |fy>| 2 , (4.D 

where p is the density of the fluid, p is the fluid pressure on the body 

surface, p^ is the pressure at infinity, and the velocity $<j> is evaluated 
on the surface of the body. 

In Figures 2a, b, we compare the pressure coefficient on a slender 
ellipsoid of revolution computed from the exact solution (see, e.g., Lamb 
[6]), the uniform perturbation solution presented in Section 2, and the panel 

method (using 300 panels). The length-to-diameter ratios for the two examples 

considered are 5 and 2.5, which correspond to values of e = 0.1 and 
e «* 0.2, respectively. We note that there is in general good agreement 
between the two approximate methods of solution and that the uniform 
perturbation method gives very good accuracy near the ends of the body, where 
non-uniform perturbation methods often fail (see, e.g., Van Dyke [7]). 
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Sorae results for the dumbbell shaped body are presented In Figures 3a, b. 

Here the perturbation method is extended only up to terms which are 
4 

0(e ), while 300 panels are still used in the panel method. The profile 
parameter b is chosen to be 3, while two different values of e are 

considered. As the figures indicate, for small values of e, the two methods 

agree well. Observe that, for these examples, as well as several other 
similar cases we considered, each perturbation solution required only about 2 
seconds of computing time, while each panel method computation required 5 
seconds (for 180 panels) or 14 seconds (for 300 panels). All computations 
were performed on a CDC Cyber 173. 

The perturbation and panel method solutions for a symmetric second— order 
Joukowski airfoil are compared for two different values of e and an angle of 
attack of 6° in Figures 5a, b. The panel method experiences some difficulty in 
resolving the flow field near the leading edge of the airfoil, while the 

uniform perturbation method produces a smooth solution. (This difficulty was 
even more evident, in other cases we considered, as the angle of attack was 
increased or the slenderness ratio decreased). In Table I, the lift 

coefficients for several different cases are calculated using these two 
methods and compared with the "exact" lift coefficient for this airfoil (see 
Karamcheti [8]). From this table, we see that the perturbation method 
produces a more accurate lift coefficient for the cases considered than does 
the panel method. 

The pressure coefficients C p on a second-order cambered Joukowski 

airfoil with two different thickness-to-chord ratios and two different angles 

of attack are presented in Figures 6a-d. The perturbation method is carried 

2 

to terms which are 0(e ), while the number of panels is 50 for each 


case. 
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The corresponding lift coefficients are summarized in Table II and compared 
with the "exact" solution for these cases. From Table II it is clear that the 
panel method produces more accurate results for thicker airfoils, while the 
uniform perturbation results are more accurate as the slenderness ratio is 
decreased or as the angle of attack is increased. This observation is 
consistent with the fact that the panel method again has difficulty in 
resolving the flow near the leading edge of the airfoil. For these examples, 
about 8 seconds of computing time were required for each perturbation 
solution, while about 5 seconds were required for each panel method 
computation. 


5. Discussion 

The examples and results displayed here have shown that the uniform 
perturbation method can be used to describe accurately some aerodynamic 
characteristics for slender bodies of revolution and for thin airfoils. When 
compared with the exact solutions available, the results of several cases 
similar to those presented here show that the perturbation solutions tend to 
become more accurate as the number of terms in the perturbation expansion is 
increased, which suggests that the perturbation solutions are convergent 
series, at least up to some maximum value of e. 

The uniform perturbation method is designed to provide an accurate 
description of the flow field over the entire surface of a smooth slender 
body, including the ends of the body. The panel method experiences some 
difficulty in a region near the leading edge of a thin airfoil (which could be 
corrected, of course, with the addition of more panels in this region or by a 
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high-order panel method). This inaccuracy in the flow description leads to a 
less accurate prediction of the lift coefficient. For thicker bodies of 
revolution or thicker airfoils, this problem is less pronounced and, in fact, 
the panel method gives more accurate results than does the (truncated) 
perturbation expansion we have used here. Thus, we find that the panel and 
uniform perturbation methods nicely complement each other, with the uniform 
perturbation solution providing a more accurate description of the flow for 

very slender bodies and the panel method being more accurate for thicker 
bodies. 

In addition to the accuracy considerations just mentioned, it should be 
emphasized that the perturbation method often requires less computer storage 
and execution time than does the panel method. The primary mathematical 
operations involved in the computation of the perturbation solution are simply 
some functional evaluations and numerical integrations. Since relatively 
little storage is required for these computations, they could be performed on 
a microcomputer, although we have made no effort to carry out our calculations 
in this manner. 

Our studies demonstrate that the uniform perturbation method can be used 
to provide low-cost, accurate solutions to flow problems involving slender 
bodies. Because of their low cost, these solutions could also be used 
effectively to provide qualitatively reliable solutions in the initial stages 
of design efforts for even thicker bodies of revolution and airfoils, with the 
panel method being used only when more accuracy is required. 
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Table I: Comparison of values for the lift coefficient for a symmetric Joukowski airfoil with different 


slenderness ratios and angles of 

attack using 

the exact, perturbation. 

and panel methods. 

Slenderness 

ratio 

e 

Exact 

C A,ex 

Perturbation 

C £,pt 

Panel 

C i,pn 

C “ c 

&,pt £,ex inn ff/ 

C — C 

£,pn *,ex „ 

p ^ lUU/o 

*,ex 

_ X 1UU/. 

x,ex 




Y - 6° 



0.01624 

0.67278 

0.67319 

0.63080 

0.06 

-6.24 

0.03248 

0.68797 

0.68961 

0.65377 

0.24 

-4.97 

0.04872 

0.70234 

0.70603 

0.66472 

0.53 

-5.36 

0.06496 

0.71588 

0.72245 

0.67684 

0.92 

-5.45 




Y = 12° 



0.01624 

1.33819 

1.33901 

1.23750 

0.06 

-7.52 

0.03248 

1.36840 

1.37167 

1.29823 

0.24 

-5.13 

0.04872 

1.39698 

1.40432 

1.32128 

0.53 

-5.42 

0.06496 

1.42391 

1.43698 

1.34485 

0.92 

-5.55 


Table II: Comparison of values for the lift coefficient for a cambered Joukowski airfoil with different 
slenderness ratios and angles of attack using the exact, perturbation, and panel methods. 


Slenderness 

ratio 

e 


0.02688 

0.05376 

0.08064 

0.10752 


0.02688 

0.05376 

0.08064 

0.10752 


0.02688 

0.05376 

0.08064 

0.10752 


Exact Perturbation Panel 

c c c C)1 -P‘ ' C *’^ x 100% ' C *> ex 

I, ex *,pt *>P n c jt >e x C £,ex 




-< 

11 

o 

o 



0.15689 

0.16101 

0.14812 

2.62 

-5.59 

0.31268 

0.32987 

0.30390 

5.50 

-2.81 

0.46641 

0.50658 

0.46610 

8.61 

-0.07 

0.61721 

0.69115 

0.63346 

11.98 

2.63 



,0 

Y = 6 



0.82841 

0.83332 

0.77878 

0.59 

-5.99 

0.99734 

1.01767 

0.95359 

2.04 

-4.39 

1.16264 

1.20984 

1.12378 

4.06 

-3.34 

1.32350 

1.40981 

1.29947 

6.52 

-1.82 



Y - 12° 



1.49084 

1.49649 

1.38619 

0.38 

-7.02 

1.67107 

1.69432 

1.58984 

1.39 

-4.86 

1.84614 

1.89984 

1.76670 

2.91 

-4.30 

2.01528 

2.11303 

1.94816 

4.85 

-3.33 


x 100% 
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Figure 2a 
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Figure 2b 
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Figure 3a 
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Figure 5a 








Figure 6b 
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Figure 6c 
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Figure Captions 

1. An axially symmetric body immersed in a uniform stream, with an indication 
of the coordinate system and the region of the singularity distribution. 


2. Pressure coefficient Cp(x) for an ellipsoid of revolution with S(x) = 
4x(l - x) and (a) e = 0.1 and (b) e = 0.2. The solid line represents 
the perturbation solution, "o" represents the exact solution, and 
" ||" represents the numerical solution. 


3. Pressure coefficient Cp(x) for a dumbbell shaped body with S(x) = 

12x(l - x)[l - 3x(l - x)] and (a) e = 0.03 and (b) e = 0.07. The solid 
line represents the perturbation solution and " |_|" represents the 
numerical solution. The profile shape of the body is indicated at the 
bottom of each figure. 

4. A two-dimensional airfoil immersed in a uniform stream, with an indication 
of the coordinate system and the endpoints (a and P) of the line 
singularity distribution. 

5. Pressure coefficient Cp(x) for a symmetric Joukowski airfoil with 
C(x) = 0 and S(x) = x(l - x) and angle of attack y " 6° for 

(a) e = 0.01624 and (b) e = 0.04872. The solid line represents the 
perturbation solution while " |~|" represents the numerical solution. 
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6. Pressure coefficient Cp(x) calculated from the perturbation solution 
(solid line) and the panel method ( |_| ) for a cambered second-order 
Joukowski airfoil with C(x) = x(l - x) , S(x) = x(l - x)^, and 
(a) y = 0°, e = 0.02688; (b) y - 0°, e = 0.08064; (c) Y = 6°, 
e = 0.02688; and (d) y - 6°, e = 0.08064. 
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